global boundaryPlot
global entropyPlot
global computeCorrC
global topological
close all
boundaryPlot = 1;
entropyPlot = 1;
computeCorrC = 0;
topological = 1;
dataMode = 1;
% dataMode = 
% 1: simulate date ; 2: flow field data ; 3: generate data ;

dirName = 'Data\1203a';
% dirName = 'G';
% dirName = 'flowFieldData\flowData1';
fileNum = 10; % in direct "G\" or "flowFieldData\flowData2", 
M = 20; N = 35;
% when topological == 1, M & N stand for nc 's range (about 15~35)
% when topological == 0, M & N stand for r 's range (about 2~5)

if entropyPlot == 1
    ncMi =plotEntropyFix(dirName,fileNum,M,N,dataMode); 
    ncMt = ncMi;
else
    ncMi = 40;
    ncMt = 40; % given ncM
end

computeCorrelationExp(dirName,fileNum,dataMode);
for ncM = ncMi : ncMt
    computeCorrelation_TheoryFix( dirName,fileNum,ncM,dataMode);
end

clear global boundaryPlot
clear global entropyPlot
clear global computeCorrC
clear global topological

% computeCorrelation_TheoryFree( dirName,fileNum,i);
